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Abstract: 

We  describe  a  methodology  for  modeling  heart  rhythms  observed  in  electrocardio¬ 
grams.  In  particular,  we  present  a  procedure  to  derive  simple  dynamic  models  that 
capture  the  cardiac  mechanisms  which  control  the  particular  timing  sequences  of 
P  and  R  waves  characteristic  of  different  arrhythmias.  Important  aspects  of  our 
models  are  their  ease  of  construction  and  conciseness.  Specifically,  these  models 
consist  of  a  structural  level,  representing  interactions  among  various  cardiac  elec¬ 
trical  events,  and  a  parameter  level,  defining  timing  statistics  of  these  events  and 
their  interactions.  The  modeling  procedure  is  a  two-step  process:  By  treating  the 
cardiac  electro-physiology  at  an  aggregate  level,  simple  network  models  of  the  wave 
generating  system  under  a  variety  of  diseased  conditions  can  be  developed.  These 
network  models  are  then  systematically  converted  to  stochastic  Petri  nets  which 
offer  a  compact  mathematical  framework  to  express  the  dynamics  and  statistical 
variability  of  the  wave  generating  mechanisms.  Models  of  several  arrhythmias  are 
included  in  order  to  illustrate  the  methodology.  One  potential  application  for  these 
models  is  in  the  development  of  automatic  classification  schemes  for  cardiac  arrhyth- 
mi»c,  as  the  models  can  be  used  as  the  basis  for  hypothesis  testing  and  parameter 
estimation  algorithms. 
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1.  Introduction 

Many  common  cardiac  arrhythmias  can  be  diagnosed  based  solely  on  the  timing 
of  various  types  of  waves  observed  in  electrocardiograms  (ECG’s).  To  develop  a 
classification  algorithm  for  cardiac  arrhythmias,  these  wave  arrival  patterns  must 
be  represented  in  a  mathematically  formal  way.  The  representation  must  be  sophis¬ 
ticated  enough  to  capture  a  wide  variety  of  wave  patterns,  yet  at  the  same  time  it 
most  be  simple  enough  to  facilitate  formulation  of  a  classification  algorithm.  In  this 
paper  we  do  not  develop  classification  algorithms  but  focus  instead  on  developing 
a  new  methodology  for  the  construction  of  concise  models  that  accurately  capture 
the  wave  timing  characteristics  of  a  wide  variety  of  cardiac  rhythms. 

The  motivation  for  our  work  comes  from  the  successes  and  limitations  of  previous 
studies.  In  particular,  in  some  previous  studies,  relatively  simple  Markov  models 
have  been  used  to  describe  patterns  of  particular  observations  of  cardiac  events. 

While  the  use  of  such  models  allows  one  to  use  powerful,  statistically  optimal  clas¬ 
sification  algorithms,  the  major  limitation  is  that  such  models  can  be  developed  and 
are  appropriate  only  for  a  limited  class  of  arrhythmias.  For  example,  three-state 
Markov  chains[l]  and  state-space  formulations [2]  have  been  used  to  model  timing 
patterns  for  the  most  prominent  ECG  feature,  the  R  wave,  and  successful  classifi¬ 
cation  algorithms  have  been  developed  based  on  them  (e.g.,  in  [2]  Kalman  filters 
were  applied  to  computation  of  the  “likelihoods”  in  a  hypothesis  testing  scheme.). 

Neither  of  these  models,  however,  can  accomodate  ECG  features  called  the  P  waves 

which  are  essential  in  the  characterization  of  many  cardiac  arrhythmias.  Moreover,  ^ 

in  some  arrhythmias  (such  as  the  Wenckebach  rhythm  described  later  in  this  pa-  (  1>0 

I _ 

per)  the  observed  wave  sequence  pattern  is  complex,  and  the  development  of  simple 
and  robust  dynamic  models  for  these  rhythms  is  not  obvious  if  one  approaches  the 
problem  purely  from  a  “black-box”  viewpoint.  n  por 

These  limitations  of  previous  models  provided  the  motivation  for  Doerschuk’s  *  at  I 
work[3]  on  “interacting  Markov  chains”  models.  In  this  work  Doerschuk  has  sue-  ,ed  ^ 

ceeded  in  developing  statistical  models  that  accurately  capture  the  timing  behavior  'tion — _ _ 

of  P  and  R  waves  together.  His  success  can  be  attributed  to  the  “physiological 

v  Distribution/ 
LjvaUablllty  Codes" 

;  'Avail  and/or 

Dlst  I  Special 


depth”  of  the  models.  While  others  tried  to  model  the  wave  patterns  directly  in 
a  black-box  fashion,  he  instead  modeled  the  physiological  events  responsible  for 
generation  of  the  ECG  waves.  As  we  will  see  in  this  paper,  the  phenomenology  of 
such  physiological  events  is  more  consistent  and  concise;  thus,  it  is  more  suitable  as 
the  foundation  of  models  than  the  wave  patterns  themselves.  The  shortcoming  of 
Doerschuk’s  models,  however,  was  their  complexity  due  to  the  nature  of  the  model 
representations  used.  In  particular,  the  desire  to  use  Markov  models,  i.e.,  the  de¬ 
sire  to  directly  construct  a  state  description,  forces  one  to  deal  simultaneously  with 
fine-level  timing  parameters  and  more  aggregate  structural  aspects  of  the  model. 
Because  of  this,  one  loses  some  of  the  conciseness  of  description  that  one  would  like 
both  for  model  construction  and  for  any  analysis  or  algorithm  design  based  on  these 
models. 

In  this  paper  we  present  a  procedure  for  modeling  the  electrical  events  char¬ 
acterizing  various  cardiac  arrhythmias  using  stochastic  Petri  nets.  We  feel  that 
stochastic  Petri  nets  are  more  suitable  for  concise  modeling  of  the  cardiac  electrical 
events  than  “interacting  Markov  chains”  and  are  capable  of  avoiding  the  complexity 
experienced  by  Doerschuk[3].  In  particular,  in  this  framework  we  can  quite  easily 
separate  and  control  the  two  significant  aspects  of  cardiac  activity  highlighted  in 
Doerschuk’s  model  -  namely  the  timing  of  events  in  different  parts  of  the  heart 
and  the  interactions  among  these  parts.  As  we  will  see,  the  interactions  specify  the 
complete  structure  of  the  Petri  net  model,  while  timing  information  affects  specific 
parameter  values  within  the  structure. 

In  Section  2  of  the  paper,  we  review  those  aspects  of  cardiac  electro-physiology 
of  importance  for  the  modeling  of  ECG  rhythms.  We  place  particular  emphasis 
on  the  mechanisms  generating  sequences  of  P  and  R  waves.  We  then  represent 
these  mechanisms  as  dynamic  systems  containing  concurrently  operating,  stochas¬ 
tic,  timing  processes.  In  Section  3,  stochastic  Petri  nets  are  introduced.  As  this 
modeling  formalism  was  expressly  developed  to  model  concurrently  operating  and 
interacting  processes,  it  is  a  natural  choice  for  modeling  cardiac  rhythm.  In  Section 
4  we  show  how  the  flow  models  developed  in  Section  2  can  be  transformed  into 
dynamical  models  of  the  type  described  in  Section  3.  In  particular,  a  systematic 
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procedure  to  generate  stochastic  Petri  net  models  of  various  cardiac  arrhythmias  is 
described.  We  illustrate  the  procedure  through  the  construction  and  simulation  of 
models  for  several  different  cardiac  arrhythmias.  Finally,  the  paper  concludes  with 
closing  discussions  in  Section  5. 
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2.  The  Cardiac  Electrical  Conduction  System 


In  this  section  we  identify  and  categorize  various  electrical  events  in  the  heart  in 
order  to  facilitate  modeling  of  cardiac  arrhythmias.  We  describe  three  important 
features  of  the  cardiac  electro-physiology:  autorhythmic  rate,  conduction  delay,  and 
refractory  period.  Based  on  these  three  we  then  define  the  basic  subunits  of  the 
electrical  conduction  system  in  the  heart  —  rhythm  and  transmission  elements.  We 
show  that  various  common  cardiac  arrhythmias  can  be  dynamically  characterized 
by  networks  of  these  rhythm  and  transmission  elements. 

2.1.  Electrical  Activity  in  the  Heart  Cells 

General  characteristics  of  the  electrically  active  cells  in  the  heart  have  been  well 
studied  and  documented.  The  following  is  a  brief  summary  of  those  characteristics 
important  in  explaining  the  mechanisms  of  many  common  cardiac  arrhythmias.  For 
more  details  see,  for  example,  [11]  and  [12]. 

The  ECG  is  a  recording  of  the  net  electrical  potential  difference  produced  when 
the  cardiac  muscle  cells  are  excited  electrically.  All  cardiac  muscle  cells  follow 
the  same  electrical  excitation  pattern:  At  rest,  they  maintain  a  certain  voltage 
(about  —90  mV)  across  the  cell  membrane.  But  if  a  sudden  change  in  the  electrical 
environment  of  the  cell  raises  the  cross- membrane  voltage  above  a  certain  threshold, 
then  the  cell  reacts  to  increase  the  voltage  spontaneously  to  a  higher  level  (about 
10  mV).  This  spontaneous,  quick  increase  in  cross-membrane  voltage  is  the  onset 
of  the  excitation.1  The  voltage  soon  falls  back  to  the  original  resting  level  (thus 
forming  an  electrical  pulse).  However,  it  takes  the  cells  a  little  longer  to  return 
to  the  same  excitable  state.  During  this  period  the  cells  cannot  be  induced  to 
increase  their  cross-membrane  voltages  spontaneously,  and  this  duration  is  called 
the  refractory  period.  The  cardiac  muscle  cells  come  out  of  the  refractory  period  in 
about  300  miliseconds  after  the  onset  of  the  excitation,  at  which  time  they  become 
ready  for  another  excitation. 

Most  commonly,  electrical  excitations  in  a  cardiac  muscle  cell  are  triggered  by 

1This  phenomenon,  observed  also  in  nerve  cells  and  skeletal  muscle  ceils,  is  called  depolarization. 


excitations  in  its  neighboring  cells.2  For  example,  if  a  section  of  a  muscle  tissue  from 
the  heart  is  electrically  stimulated,  the  excitation  spreads  quickly  into  the  entire 
tissue.  In  such  inter-cellular  transmission  of  electrical  excitation,  the  refractory 
period  is  important  because  it  prevents  back-flows  and  reflections.  For  example, 
suppose  that  the  excitation  of  cell  A  triggers  excitation  in  cell  B,  a  neighbor.  Since 
cell  A  will  be  in  refractory  period,  cell  B  will  not  be  able  to  transmit  the  same 
excitation  back  to  cell  A.  But  suppose  that  cell  B  has  a  neighbor  cell  C  which  is 
not  a  neighbor  of  cell  A  and  is  ready  to  be  excited.  Then  cell  B  can  transmit  the 
excitation  to  cell  C.  The  refractory  period,  thus,  gives  directionality  in  inter-cellular 
transmission  of  electrical  excitation. 

As  described  above,  most  cardiac  muscle  cells  are  induced  into  activity  by  re¬ 
ceiving  excitation  from  neighboring  cells.  However,  many  cardiac  muscle  cells  also 
have  a  mechanism  to  generate  periodic  excitations  by  themselves.  In  fact,  the  cross¬ 
membrane  voltage  in  these  cells  at  rest  rises  slowly  towards  the  threshold  for  the 
spontaneous,  quick  increase  that  characterizes  the  excitation  of  the  cells.  This  prop¬ 
erty  is  called  autorhythmicity.  The  rate  at  which  the  resting  voltage  reaches  the 
excitation  threshold  differs  from  cell  to  cell.  The  heart  has  a  collection  of  cells  which 
have  the  fastest  of  such  rates  and  specialize  in  the  generation  of  periodic  excitations. 
In  a  normal  heart,  the  cells  in  a  region,  called  the  SA  node  on  the  wall  of  the  right 
atria  are  the  source  of  these  excitations.  The  periodic  electrical  activity  of  the  SA 
node  is  eventually  transmitted  throughout  the  heart,  and  the  muscle  cells  of  the 
atria  and  ventricles  contract,  just  as  all  skeletal  muscle  cells  do,  upon  reception  of 
the  electrical  excitation.  The  rate  at  which  the  SA  node  excites  is,  thus,  observed 
as  the  heart  rate  in  a  normal  heart. 

2.2.  Normal  Heart  Activity  and  the  Electrocardiogram 

The  heart  consists  of  two  pairs  of  chambers:  the  atria  and  vcntnelet.  The  two  pairs^ 
contract  alternately  in  such  a  fashion  that  the  atria  pump  blood  from  the  veins 
into  the  ventricles  which  then  pump  it  to  the  lung  and  the  rest  of  the  body  At 

3This  is  accomplished  through  the  gap  junction*,  croee-csllul ai  chstmti*  for  >ntr*  <s lluUi  plasm.. 
Increase  in  electrical  potential  inside  a  cell  causes  the  potentials  instd*  its  nsisM.-.r  i.iu  i<  r i»» 
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contraction,  an  electrical  excitation  is  transmitted  quickly  throughout  the  muscle 
cells  of  the  chamber  pair.  This  nearly-simultaneous  excitation  of  many  cells  is 
observed  as  a  peak  on  the  ECG  recording.  Because  the  ventricles  have  much  larger 
muscle  mass  than  the  atria,  they  produce  a  conspicuous  peak  called  the  R  wave3, 
while  the  atria  produce  a  small  (and  often  difficult  to  detect)  peak  called  the  P  wave. 
The  timing  of  chamber  contractions  is  orchestrated  by  a  network  of  muscle  cells 
specialized  in  generation  and  transmission  of  electrical  excitations.  The  electrical 
activities  in  a  normal  heart  are  as  follows:  The  SA  node  generates  excitations 
periodically  due  to  its  autorhythmicity4.  These  excitations  spread  throughout  the 
atria,  inducing  contractions.  They  also  reach  the  AV  node  on  the  wall  separating 
the  atria  from  the  ventricles.  The  AV  node  is  the  only  normal  electrical  channel 
between  the  two  pairs  of  heart  chambers,  and  its  excitation  is  normally  transmitted 
to  the  muscle  cells  of  the  ventricles  through  specialized  conducting  muscle  cells 
called  the  Purkinje  fibers.  Unlike  intra-muscluar  conduction,  the  transmission  of 
electrical  excitation  through  the  AV  node  is  slow.  In  effect,  this  delays  contraction 
of  the  ventricles  from  that  of  the  atria,  so  that  the  blood  can  flow  in  the  correct 
direction  from  the  atria  to  ventricles.  The  delay  through  the  AV  node,  thus,  plays 
an  important  role  in  the  control  of  chamber  contraction. 

In  rhythm  analysis  of  the  ECG,  we  try  to  make  an  assessment  on  how  well 
this  electrical  conduction  system  is  controlling  chamber  contraction.  The  intervals 
between  P  and  R  waves  are  studied  for  this  purpose.  In  many  diseases,  parts 
of  the  cardiac  electrical  conduction  system  are  affected,  and  we  observe  distinct, 
anomalous  patterns  in  P,R  wave  sequences.  These  known  patterns  form  a  very 
important  diagnostic  basis  for  cardiac  arrhythmias. 

2.3.  Control  of  the  Cardiac  Electrical  Conduction  System 

The  electrical  conduction  system  of  the  heart  has  two  main  functions:  generation 
of  electrical  excitations  and  distribution  of  these  throughout  the  heart.  Generally 

3Strictly  speaking  this  should  be  the  QRS  complex,  a  sequence  of  three  waves.  But  for  simplicity, 
we  use  the  term  *R  wave*  to  refer  to  the  ECG  feature  corresponding  to  the  contraction  of  the 
ventricles. 

4The  autorhythmic  rate  of  the  SA  node  is  under  autonomic  nervous  control. 
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speaking,  these  functions  are  regulated  by  three  timing  parameters:  autorhythmic 
rate,  conduction  delay,  and  refractory  period.  This  subsection  describes  how  these 
timing  parameters  affect  the  functions  of  the  system.  1)  Autorhythmic  rate.  Nor¬ 
mally,  the  SA  node  has  the  fastest  autorhythmic  rate  in  the  whole  system,  and  the 
heart  activities  are  paced  at  this  rate.  However,  non-SA  nodal  tissue  can  sometimes 
attain  an  autorhythmic  rate  comparable  to  that  of  the  SA  node.  In  such  a  system, 
two  or  more  rhythm  sources  compete  to  pace  the  heart,  resulting  in  an  anomalous 
and  possibly  chaotic  contraction  pattern.  This  phenomenon  can  occur  either  when 
the  rate  at  the  SA  node  slows  down  abnormally  or,  more  commonly,  when  the  rate 
of  non-SA  nodal  tissue  speeds  up  drastically.  2) Conduction  delay.  As  mentioned 
before,  the  conduction  speed  through  the  AV  node  has  an  important  function  in 
delaying  ventricular  contractions  after  the  preceding  atrial  contractions.  This  is  an 
example  of  the  direct  influence  that  conduction  delays  generally  have  over  cardiac 
control.  3) Refractory  period.  An  abnormally  long  refractory  period  can  period¬ 
ically  block  excitations.  For  example,  suppose  that  the  duration  from  the  onset 
of  excitation  to  the  end  of  the  refractory  period  at  the  AV  node  is  slightly  longer 
than  the  period  of  autorhythmic  excitations  at  the  SA  node.  The  first  excitation 
from  the  SA  node  can  successfully  excite  the  AV  node.  The  second  excitation, 
however,  arrives  at  the  AV  node  during  the  refractory  period,  and  it  is  blocked 
from  propagating  further.  The  resulting  P,R  sequence  in  this  hypothetical  case  is  a 
series  of  P  waves  alternatingly  followed  by  R  waves.  In  reality,  blockage  of  excita¬ 
tion  by  the  refractory  periods  is  not  absolute.  In  general,  the  refractory  period  is 
loosely  divided  into  two  parts.  The  early  part  is  where  excitation  is  blocked  with 
certainty  and  is  called  the  absolute  refractory  period.  In  the  latter  part,  the  tissue 
goes  through  a  period  where  it  becomes  increasingly  prone  to  be  excited,  and  this 
period  is  called  the  relative  refractory  period.  Forcing  tissues  to  excite  during  their 
relative  refractory  period  often  alters  their  autorhythmic  and  trasmission  charac¬ 
teristics  momentally.  For  example,  an  excitation  arriving  at  the  AV  node  during  its 
relative  refractory  period  will  take  longer  than  usual  to  propagate  through  the  node. 
Also,  if  a  pacing  tissue  like  the  SA  node  for  some  reason  receives  an  excitation  dur¬ 
ing  a  relative  refractory  period,  its  autorhythmic  generation  of  the  next  excitation 
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pulse  is  momentally  delayed.  Such  brief  interference  with  natural  autorhythmicity 
is  referred  to  as  either  the  resetting  or  stunning  phenomenon.  In  resetting,  the  time 
interval  between  the  reception  of  the  external  excitation  and  the  generation  of  the 
next  autorhythmic  excitation  is  roughly  equal  to  the  period  of  the  autorhythmicity. 
In  stunning,  on  the  other  hand,  this  time  interval  is  significantly  longer  than  the 
period  of  the  autorhythmicity. 

2.4.  The  Rhythm  and  Transmission  Elements 

Since  the  two  main  functions  of  the  cardiac  electrical  conduction  system  are  genera¬ 
tion  and  distrubution  of  electrical  excitations,  the  dynamics  of  the  system  can  be  de¬ 
scribed  by  a  network  of  two  types  of  elements.  One  of  these  is  the  rhythm  element, 
which  generates  electrical  excitations,  and  the  other  is  the  transmission  element, 
which  distributs  excitation  from  one  section  of  the  system  to  another.  The  control¬ 
ling  mechanisms  for  macroscopic  flows  of  electrical  signals  (excitations)  in  the  heart 
can  be  represented  by  a  network  of  rhythm  and  transmission  elements.  In  the  next 
subsection,  we  represent  the  dynamics  of  several  cardiac  arrhythmias  with  such 
networks.  In  this  subsection  the  two  elements  are  described. 

We  represent  rhythm  elements  diagramatically  as  triangles  and  transmission 
elements  as  rectangles  (Fig.  1).  The  rhythm  elements  are  used  to  describe  the 
autorhythmic  properties  of  cardiac  muscle  tissues;  their  primary  function  is  periodic 
generation  of  electrical  signals.  Associated  with  a  rhythm  element  are  input  and 
output  terminals  as  well  as  several  parametric  variables  that  control  the  intervals  of 
signal  generation.  The  purposes  of  these  terminals  and  parameters  are  as  follows: 
(i)  output.  The  signal  generated  by  the  element  is  sent  to  neighboring  elements 
through  the  output  terminal,  (ii)  input.  The  autorhythmic  ceils  can  be  stimulated 
by  excitations  in  neighboring  cells.  The  rhythm  element  receives  such  external 
stimulation  through  the  input  terminal,  (iii)  variables.  The  fundamental  parameter 
of  a  rhythm  element  is  the  period  at  which  it  generates  signals.  However,  this 
basic  period  can  be  altered  especially  when  an  external  stimulation  arrives  through 
the  input  terminal.  Thus,  besides  the  basic  period,  variables  that  quantitatively 
characterize  the  external  influence  on  the  function  of  the  element  are  needed.  For 


l 


>.Y>. 


example,  one  such  variable  is  the  refractory  period.  Others,  associated  with  such 
phenomena  as  resetting  and  stunning,  will  be  described  in  Section  4. 

The  transmission  elements  are  used  to  describe  the  delay  of  signals  traveling 
through  cardiac  muscle  tissues.  It  is  basically  a  bi-dirctional  channel  characterized 
by  two  pairs  of  input  and  output  terminals.  Signals  can  be  transmitted  through  the 
element  in  either  direction,  but  when  two  opposing  signals  meet  in  the  element  they 
annihilate  each  other.  The  variables  associated  with  the  transmission  elements  are 
the  transmission  delays  for  both  directions,  refractory  periods  which  characterize 
the  excitability  of  the  two  input  terminals,  and  other  parameteric  variables  repre¬ 
senting  the  factors  that  may  influence  the  durations  of  transmission  delays  and/or 
refractory  periods.  (We  will  show  an  example  of  such  factors  later  when  we  describe 
the  Wenckebach  phenomenon.) 

The  rhythm  and  transmission  elements  are  used  quite  flexiblly,  and  in  particular, 
the  numbers  and  types  of  the  terminals  and  variables  assigned  to  an  element  are 
adjustable.  For  example,  a  uni-directional  transmission  element  would  graphically 
be  presented  as  a  rectangle  just  like  a  bi-directional  transmission  element(Fig.  lb) 
but  without  the  second  pair  of  input  and  output  terminals.  Such  variations  among 
rhythm  and  transmission  elements  are  presented  in  more  details  in  Section  4,  when 
we  discuss  the  implementation  of  these  elements  with  stochastic  Petri  nets. 

2.5.  Signal  Flow  Block  Diagrams 

The  rhythm  and  transmission  elements  described  in  the  preceding  subsection  allow 
us  to  model  the  timings  of  P  and  R  waves  and  aberrancies  of  different  rhythms  at  a 
relatively  aggregate  level.  While  it  is  certainly  possible  to  use  the  modeling  method¬ 
ology  developed  in  this  paper  to  describe  cardiac  activity  at  a  more  detailed  level 
(by  partitioning  the  heart  into  a  larger  number  of  interacting  rhythmic  and  conduc¬ 
tive  units,  each  representing  a  smaller  portion  of  cardiac  electrical  pathways),  the 
explicit  use  to  which  we  put  this  methodology  here  is  at  the  other  extreme.  Specif¬ 
ically,  by  highlighting  the  mechanism  causing  and  driving  particular  arrhythmias, 
we  want  to  obtain  the  simplest  possible  models  capturing  characteristics  of  the  cor¬ 
responding  P,R  sequences.  Such  “minimal  representations”  should  ultimately  be 
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of  most  value  as  the  basis  for  robust  signal  processing  and  automated  diagnosis  of 
cardiac  arrhythmias. 

To  illustrate  this  philosophy,  consider  the  modeling  of  a  perfectly  normal  heart¬ 
beat  sequence.  We  can  divide  the  cardiac  condution  system  into  five  stages  based 
on  their  structural  and  functional  differences  —  the  SA  node,  intra-atrial  conduc¬ 
tive  paths,  AV  node,  Purkinje  fiber  conductive  paths,  and  ventricles.  Each  tissue 
block  can  in  fact  excite  autorhythmically,  and  all  except  the  SA  node  and  ventricles 
(which  are  the  two  ends  of  this  electrical  system)  can  conduct  bi-directionally.  As 
far  as  the  modeling  of  a  normal  P,R  sequence  is  concerned,  however,  bi-directional 
conduction  is  an  unnecessary  physiological  detail  since  excitations  conducting  in 
the  retrograde  direction  (i.e.,  direction  towards  the  SA  node,  opposite  to  the  nor¬ 
mal  direction  of  the  flow  of  excitation)  are  never  observed  in  such  a  sequence. 
Thus,  a  model  constructed  with  uni-directional  transmission  elements  instead  of 
bi-directional  elements  is  simpler  yet  phenomenologically  just  as  acurate.  Such  a 
signal  flow  block  diagram  model  is  shown  in  Figure  2a.  The  rhythm  element  rep¬ 
resenting  the  SA  node  does  not  have  an  input  terminal  because  of  the  absence  of 
retrograde-conducting  excitation.  The  autorhythmic  excitations  of  the  SA  node 
activate  the  atria  and  generate  the  P  waves.  The  letter  “P”  by  the  output  terminal 
represents  the  generation  of  the  P  wave.  The  excitation  also  bifurcates  to  activate 


the  trasmission  element  and  reset  the  rhythm  element  of  the  intra-atrial  pathways. 
The  excitation  produced  by  the  intra-atrial  pathways  is  the  result  of  either  the 
transmission  of  the  excitation  originated  in  the  SA  node  or  its  own  autorhythmic 
activation.  This  is  represented  by  the  convergence  of  the  outputs  of  the  transmission 
and  rhythm  elements.  The  AV  node  and  Purkinje  fibers  are  reprsented  by  a  par¬ 
allel  pair  of  transmission  and  rhythm  elements,  just  like  the  intra-atrial  pathways. 
The  only  differences  among  these  three  parts  of  the  block  diagram  are  the  values  of 
the  parameters  assigned  to  the  respective  transmission  and  rhythm  elements  (i.e., 
conduction  delays,  autorhythmic  intervals,  and  absolute  refractory  periods).  The 
output  of  the  Purkinje  fibers  is  directed  to  the  ventricles  whose  activation  produce 
the  R  waves.  The  letter  “R”  by  the  input  of  the  ventricles  represents  the  production 
of  the  waves.  This  block  diagram  model  implies  that  the  autorhythmic  period  of  the 


ventricles  is  larger  than  those  of  other  parts  of  the  heart  (in  particular  the  SA  node) 
so  that  the  ventricles  are  “reset”  frequently  enough  not  to  excite  autorhythmically. 
In  the  model,  therefore,  the  rhythm  element  representing  the  ventricles  does  not 
have  an  output  terminal.  This  is  consistent  with  the  assumption  that  retrograde 
conduction  is  absent. 

This  block  diagram  model  of  the  normal  heart  can  be  simplified  further  as 
follows:  First,  since  in  a  normal  cardiac  sequence  the  autorhythmic  rate  in  the  SA 
node  is  the  fastest  and  drives  the  entire  system,  the  rhythm  elements  in  the  rest  of 
the  system  are  always  “reset”  before  they  can  generate  an  autorhythmic  excitation. 
We  can,  therefore,  delete  the  three  rhythm  elements  used  to  model  the  autorhythmic 
properties  of  the  intra-atrial  pathways,  AV  node,  and  Purkinje  fibers.  After  these 
three  rhythm  elements  are  removed,  the  three  remaining  stages  that  separate  the 
two  wave  generators  collectively  form  only  a  series  of  three  transmission  elements. 
These  transmission  elements  are,  thus,  combined  into  one  aggregate  element,  which 
we  call  the  “AV  node”  for  simplicity.  The  simplified  model  of  the  normal  heart 
is  shown  in  Figure  2b,  which  has  only  two  rhythm  elements  and  one  transmission 
element. 

2.6.  Modeling  of  Common  Arrhythmic  ECG  Patterns 

In  this  subsection,  we  represent  the  dynamics  of  several  common  cardiac  arrhyth¬ 
mias  using  signal  flow  block  diagrams  made  up  of  rhythm  and  transmission  elements. 
Many  of  the  common  arrhythmias  are  included  in  the  four  categories  described  in 
this  subsection,  and  Figures  3a  to  8a  show  signal  flow  models  of  several  examples 
from  each  category.  Also,  to  illustrate  the  wave  patterns  for  these  examples  of 
arrhythmias,  Figure  3b  to  8b  present  the  P,R  wave  sequences  obtained  from  simu¬ 
lations  of  the  models  (See  Section  4  for  details).  P  waves  are  represented  by  short 
vertical  lines  and  R  waves  by  long  vertical  lines,  and  abnormal  waves  are  repre¬ 
sented  by  lines  with  small  squares  at  their  tips.  The  interval  of  time  ticks  is  one 
second.  See  [10]  for  more  detailed  descriptions  of  cardiac  arrhythmias. 

1.  Extrasystole  -  ectopic  beat 

(i)phenomtnology.  In  this  condition,  non-SA  nodal  tissue  becomes,  either  contin- 


uously  or  sporadically,  a  pacemaker  for  the  heart.  Such  abnormal  pacemakers  are 
called  ectopic  rhythm  sources.  Ectopic  beats  may  originate  in  the  atrial  wall  outside 
the  SA  node  causing  untimely  atrial  contractions  and  abnormally  shaped  P  waves. 
This  class  of  cardiac  arrhythmia  is  called  Atrial  Premature  Beat,  or  APB.  Ectopic 
beats  may  also  arise  in  the  AV  node,  in  the  Purkinje  fibers,  or  in  the  ventricular 
musculature,  and  they  tend  to  cause  Ventricular  Premature  Beat,  or  VPB.  These 
ectopic  beats,  in  addition,  may  initiate  excitations  that  flow  backward  into  the  atria 
(i.e.,  retrograde  conduction).  The  resulting  atrial  electrical  activity  is  referred  to 
as  the  retrograde  P  wave.  As  mentioned  before,  the  SA  node  can  be  “reset”  or 
“stunned”  by  retrograde  excitations,  (ii)modeling.  Ectopic  beats  are  caused  by 
abnormally  fast  autorhythmic  rates  at  non-SA  nodal  tissues.  Rhythm  elements  are 
used  to  represent  such  ectopic  sources.  (Hi) examples.  Figures  3a  and  4a  show  ex¬ 
amples  of  signal  flow  block  diagrams  of  APB  and  VPB,  respectively.  Both  models 
consist  of  three  stages  -  the  atria,  AV  nodes  and  ventricles  -  just  like  the  normal 
heart  model  in  the  previous  subsection.  In  the  APB  model(Fig.  3a)  the  atrial  stage 
has  two  rhythm  elements,  one  for  the  normal  SA  nodal  beats  and  the  other  for  the 
atrial  ectopic  beats.  The  output  of  each  rhythm  element  resets  the  other  rhythm 
element  and  activates  the  transmission  element  representing  the  AV  node.  The  rest 
of  the  model  is  identical  to  the  normal  heart  described  before.  The  P  waves  pro¬ 
duced  by  atrial  ectopic  beats  usually  have  abnormal  shapes;  in  the  model  a  “P”  with 
an  overbar  denotes  that  the  output  of  the  rhythm  element  representing  the  ectopic 
source  generates  such  abnormal  P  waves.  In  the  VPB  model  (Fig.  4a)  the  ectopic 
source  is  generally  thought  to  be  in  the  ventricular  tissue.  Thus,  the  rhythm  ele¬ 
ment  representing  the  ventricles  has  an  output  terminal  whose  activation  produces 
a  premature,  abnormally  shaped  R  wave  (denoted  by  an  “R”  with  an  overbar  in  the 
figure).  Since  the  excitation  generated  by  the  ventricles  conducts  in  the  retrograde 
direction  towards  the  atria,  a  bi-dirctional  transmission  element  is  used  to  represent 
the  AV  node.  Retrograde  activation  of  the  atria  causes  a  retrograde  P  wave  (“P” 
with  an  overbar)  and  resets  the  SA  nodal  autorhythmic  source. 

2.  Extrasystole  -  coupled  beat 

(i)phenomenology.  This  category  of  cardiac  arrhythmias  also  deals  with  prema- 
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ture  waves.  Although  the  occurence  of  premature  waves  in  the  previous  category 
seem  independent  from  the  timing  of  the  normal  waves  (thus  they  are  thought 
to  be  caused  by  ectopic  rhythm  sources  acting  independently  from  the  SA  nodal 
source),  the  premature  waves  in  this  category  show  some  correlation  with  the  SA 
nodal  beats.  Specifically,  these  premature  waves  are  synchronized  with  the  normal 
beats,  and  the  intervals  between  premature  waves  and  preceding  normal  waves  are 
fairly  constant.  Such  intervals  are  called  the  coupling  intervals,  (ii)modtling.  Al¬ 
though  there  are  several  physiological  or  anatomical  explanations  for  the  origin  of 
the  coupling  intervals,  a  convenient  way  to  describe  the  phenomenon  is  to  use  a  hy¬ 
pothetical  electrical  conduction  channel  called  the  reentrant  pathway.  Conceptually, 
the  reentrant  pathway  resides  within  the  muscular  wall  of  a  heart  chamber.  It  has 
a  conduction  delay  whose  value  is  equal  to  that  of  the  coupling  interval.  It  receives 
an  excitation  when  the  chamber  is  excited  to  contract,  delays  the  excitation  for  the 
amount  of  the  time  specified  by  its  conduction  delay  parameter,  and  then  returns 
the  excitation  back  to  the  chamber.  This  leads  to  a  second,  abnormal  contraction 
of  the  chamber  (and  generation  of  an  associated  ECG  wave)  following  the  normal 
contraction  arising  from  the  direct  excitation  pathway.  In  the  block  diagram,  the 
coupling  intervals  are  described  by  uni-directional  trasmission  elements  represent¬ 
ing  the  reentrant  pathways,  (iii)example.  Figure  5a  shows  the  block  diagram  for  a 
condition  called  bigeminy,  where  a  premature  ventricular  contraction  occurs  every 
other  beat.  The  first  two  stages  of  the  model  -  the  SA  and  AV  nodes  -  are  the 
same  as  the  normal  heart.  The  output  of  the  AV  node  activates  both  the  ventricles 
(producing  the  normal  R  wave)  and  the  reentrant  pathway.  The  reentrant  path¬ 
way  delays  the  excitation  and  activates  the  ventricles  (producing  the  abnormal  R 
wave).  The  interval  between  the  normal  and  succeeding  abnormal  R  waves  is  the 
coupling  interval  in  this  case,  and  it  is  modelled  by  the  delay  parameter  of  the  trans¬ 
mission  element  representing  the  reentrant  pathway.  The  normal  conduction  wave 
(initiated  by  the  SA  node)  immediately  following  the  premature  R  wave  arrives  at 
the  ventricle  during  its  absolute  refractory  period.  Thus,  the  normal  excitation  is 
blocked,  and  the  corresponding  R  wave  is  not  produced.  This  is  observed  as  the 
relatively  long  interval  between  a  premature  R  wave  and  the  succeeding  normal  R 


wave(Fig.  5b). 

3.  AV  Conduction  Block 

(i)phenomenology.  This  category  includes  arrhythmias  caused  by  abnormalities  in 
conduction  between  the  atria  and  ventricles.  In  a  complete  blockage,  called  third 
degree  AV  block,  the  AV  node  is  entirely  unable  to  conduct  excitations.  Thus,  the 
contractions  of  the  ventricles  must  be  paced  by  the  autorhythmic ity  of  the  AV  node, 
Purkinje  fibers,  or  ventricular  musculature  itself.8  The  resulting  P,R  wave  sequence 
displays  a  case  of  A  V  dissociation,  a  phenomenon  in  which  the  rhythm  of  the  R  waves 
is  independent  from  that  of  the  P  waves.  In  second  degree  AV  block,  not  all  the  atrial 
excitations  are  blocked  by  the  AV  node.  When  only  two  out  of  three  P  waves  are 
followed  by  R  waves,  the  condition  is  referred  to  as  3:2  block.  Other  ratios  commonly 
seen  are  2:1,  3:1,  4:1,  and  so  on.  A  P,R  sequence  in  this  subcategory  may  also  exhibit 
a  Wenckebach  phenomenon  described  as  follows:  The  SA  node  generates  excitations 
at  a  constant  rate,  but  the  P-R  interval  grows  progressively  longer  during  several 
beats  until  there  is  a  P  wave  not  followed  by  an  R  wave  (i.e.,  a  block  has  occurred). 
The  next  P  wave  is  followed  by  an  R  wave  with  a  normal,  short  P-R  interval;  then, 
the  interval  again  grows  progressively  longer  over  the  next  several  beats  as  the 
pattern  is  repeated.  Figure  7b  shows  an  example  of  a  P,R  sequence  displaying  the 
Wenckebach  phenomenon.  Finally,  in  first  degree  AV  block,  no  R  wave  is  missed 
after  any  P  wave,  but  the  P-R  interval  is  abnormally  prolonged,  (ii)modeling.  AV 
conduction  blocks  are  caused  by  various  disease  conditions  inside  the  AV  node. 
Since  all  the  abnormalities  occur  within  the  transmission  element  representing  the 
AV  node,  the  signal  flow  block  diagram  is  essentially  the  same  as  that  of  the  normal 
heart  (Fig.  2b).  The  transmission  element  (the  AV  node),  however,  tends  to  have 
a  considerably  more  complex  role  than  in  the  normal  P,R  sequence.  (Modeling 
of  the  third  degree  AV  block  is  an  exception  and  is  trivially  simple;  it  can  be 
accomplished  by  removing  the  transmission  element  from  the  normal  heart  model.) 
(Hi) examples.  Figure  6a  shows  a  model  of  second  degree  AV  block.  Note  that  the 
trasmission  element  has  an  extra  parameter  called  “probability  of  conduction.”  In 

6This  is  an  example  of  autorhythmic  beats  from  a  non-SA  nodal  tissue  preventing  a  complete 
failure  of  the  heart,  and  these  beats  are  called  the  escape  beats. 
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3:2  block,  for  example,  this  probability  is  set  to  |.  Figure  7a  shows  a  model  featuring 
the  Wenckebach  phenomenon.  A  completely  different  transmission  element  must  be 
used  to  describe  the  dynamics  of  the  AV  node.  In  Wenckebach,  events  in  the  current 
beat  are  dependent  on  the  events  in  the  previous  beats.  For  example,  every  P-R 
interval  must  be  longer  than  the  previous  interval  unless  the  R  wave  is  missing  from 
the  previous  beat;  if  the  R  wave  is  missed  the  P-R  interval  must  become  short  again. 
The  transmission  element,  therefore,  must  have  a  “memory”  of  its  actions  in  the 
previous  beats  to  determine  how  it  behaves  in  the  current  beat.  This  requires  us  to 
describe  the  transmission  element  itself  as  a  dynamic  system.  Detailed  specification 
of  the  behavior  of  such  a  transmission  element  is  discussed  in  Section  4  in  which  we 
describe  the  implementation  of  the  dynamics  of  various  rhythm  and  transmission 
elements. 

4.  Abnormal  conduction  path 

(i)phenomenon.  In  a  normal  heart,  electrical  insulation  exists  between  the  atria 
and  ventricles,  and  the  AV  node  is  the  only  electrical  channel  between  the  pairs  of 
chambers.  Some  abnormal  hearts,  however,  have  conduction  pathways  that  bypass 
the  AV  node.  For  example,  in  a  condition  called  the  Wolff-Parkinson- White  syn¬ 
drome  the  ECG  displays  an  abnormally  early  onset  of  the  ventricular  activity.  It 
also  exhibits  broadened  R  waves  because  each  beat  reaches  the  ventricles  through 
the  abnormal  and  normal  (AV  nodal)  routes  and  excites  different  parts  of  the  cham¬ 
bers  at  slightly  different  times,  (ii)modeling.  An  abnormal  conduction  channel  is 
represented  by  a  transmission  element,  (iii)txamph.  Figure  8a  shows  a  model  of  the 
Wolff-Parkinson- White  syndrome.  The  ventricles  are  divided  into  two  parts.  One 
of  these  is  excited  normally  and  produces  R  waves.  The  other  part  is  excited  via 
the  abnormal  conduction  channel6  and  produces  characteristic  delta  waves7.  The 
output  of  the  SA  node  simultaneously  activates  the  normal  and  abnormal  channels. 
The  conduction  delay  through  the  abnormal  channel  is  slightly  shorter  than  that 
through  the  normal,  AV  nodal  channel.  To  increase  the  anatomical  accuracy  of 
the  model,  the  two  rhythm  elements  representing  the  two  parts  of  the  ventricles 

•called  the  bundle  of  Kent 

7The  R  and  delta  waves  are  actually  parts  of  a  single  peak.  The  delta  wave  is  identifiable  as  an 
abnormally  smeared  leading  edge  of  the  conspicuous  R  wave  peak. 
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may  be  connected  with  a  bi-directional  transmission  element.  But  since  the  two 
parts  of  the  ventricles  excite  at  about  the  same  time,  such  a  transmission  element 
would  host  only  collision  and  annihilation  of  two  excitations  flowing  in  the  opposite 
directions.  This  activity  does  not  influence  the  P,R  wave  sequence  produced  by  the 
model.  The  bi-directional  transmission  element  is,  therefore,  omitted  in  favor  of 
model  simplicity. 

In  the  preceeding  description  of  the  cardiac  electrical  conduction  system,  we 
have  identified  several  of  its  dynamic  properties  which  enable  us  to  characterize 
the  P,R  wave  sequences  of  the  cardiac  arrhythmias.  In  the  following  sections,  we 
describe  such  dynamics  in  a  mathematically  formal  way  by  implementing  the  sig¬ 
nal  flow  block  diagrams  developed  in  this  section  using  stochastic  Petri  nets.  This 
process  is  similar  to  translating  flow-charts  into  computer  codes  using  a  particular 
programming  language.  It  is  desireable  that  the  programming  language  has  features 
that  allow  a  straightforward  representation  of  the  sequential  behavior  characterized 
by  the  flow-chart.  In  this  context,  two  important  features  of  the  signal  flow  block 
diagrams  for  cardiac  arrhythmias  should  be  noted:  One  is  that  the  system  dynamics 
are  regulated  by  timing  parameters,  and  the  other  is  the  concurrency  of  the  dynam¬ 
ics  (e.g.,  the  SA  node  is  in  its  refractory  period  while  the  AV  node  is  conducting 
excitation).  In  other  words,  using  the  block  diagrams  introduced  in  this  section 
we  have  characterized  the  mechanisms  of  various  cardiac  arrhythmias  as  concurrent 
timing  processes.  Stochastic  Petri  nets,  which  are  described  in  the  next  section, 
have  been  developed  to  model  dynamics  of  such  concurrent  timing  processes.  They, 
therefore,  offer  a  natural  framework  for  the  translation  process  discussed  above. 
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3.  Stochastic  Petri  Nets 

In  the  preceding  section  we  have  seen  that  the  electrical  activities  of  various  parts 
of  the  heart  are  governed  by  local  timing  parameters  such  as  autorhythmic  periods, 
conduction  delays,  and  refractory  periods.  These  activities  interact  with  each  other 
by  passing  and  receiving  excitations,  collectively  defining  the  electrical  phenomena 
observed  as  ECG’s.  Thus,  behaviors  of  the  cardiac  electrical  conduction  system  can 
be  considered  as  a  collection  of  concurrently  operating  timing  processes.  Stochastic 
Petri  nets  offer  a  general  and  flexible  format  to  express  the  activities  of  concurrent 
timing  processes.  Historically,  stochastic  Petri  nets  were  developed  as  the  result 
of  an  extension  of  the  original  Petri  net  theory  which  features  simpler  execution 
policies  (i.e.,  dynamics).  In  this  section  we  first  review  the  main  features  of  the 
original  Petri  rets  before  discussing  those  of  stochastic  Petri  nets.  The  section  closes 
with  a  discussion  of  features  of  stochastic  Petri  nets  desirable  for  the  modeling  of 
the  cardiac  electrical  conduction  system. 

3.1.  Petri  Nets 

Petri  nets  are  abstract  models  of  information  flow,  and  they  are  particularly  useful 
in  describing  and  analyzing  the  flow  of  information  and  control  in  systems  which 
exhibit  asynchronous  and  concurrent  activities.  In  this  subsection,  we  try  to  high¬ 
light  some  fundamental  properties  of  Petri  nets.  Refer  to  Peterson[4,5]  for  more 
detailed  discussions  of  them. 

Structure.  Petri  nets  are  commonly  represented  pictorally  as  directred  graphs. 
Figure  9  shows  an  example  of  a  Petri  net.  There  are  four  structual  components  — 
two  types  of  nodes,  called  places  and  transitions,  and  two  types  of  directed  arcs, 
input  and  output  arcs.  A  place  is  represented  by  a  circle  (“pO”  to  “p6”  in  Fig.  9) 
and  a  transition  by  a  bar  (“tO”  to  “t5”).  All  the  directed  arcs  in  a  Petri  net  connect 
a  node  of  one  type  with  a  node  of  the  other  type;  no  arc  connects  two  nodes  of 
the  same  type.  An  arc  going  from  a  place  to  a  transition  is  referred  to  as  either  an 
output  arc  of  the  place  or  an  input  arc  of  the  transition.  Similarly,  an  arc  going 
from  a  transition  to  a  place  is  called  either  an  output  arc  of  the  transition  or  an 
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input  axe  of  the  place. 

Dynamics.  The  dynamics  of  Petri  nets  are  represented  by  the  positions  and  move¬ 
ments  of  markers  called  tokens.  Tokens  are  pictorally  represented  by  dots  inside 
places.  In  the  example  (Fig.  9),  each  of  the  places  pO,  p2,  and  p6  has  a  single  token. 
Tokens  can  be  moved  to  other  places  along  directed  arcs  by  firing  transitions  on 
the  arcs.  A  transition  can  fire  only  if  it  is  enabled.  A  transition  becomes  enabled 
when  all  its  input  arcs  are  connected  to  places  possessing  tokens.  In  the  example, 
tl  and  t2  are  enabled,  and  t3  is  not  enabled  because  p3  does  not  have  a  token. 
Decision  Rules.  Note  that  in  the  preceding  example  tl  and  t2  cannot  fire  at 
the  same  time  because  firing  one  of  them  takes  the  token  out  of  pO  and  disables 
the  other.  Those  transitions  competing  for  the  same  token(s)  are  described  to  be 
“in  conflict."  Policies,  called  the  decision  rules,  to  resolve  such  conflicts  are  thus 
specified.  One  implementation  of  the  decision  rule  is  a  probabilistic  policy.  For 
example,  in  Figure  9  we  can  make  t2  fire  with  probability  of  0.7  whenever  tl  and 
t2  are  in  conflict.  Note  that  when  this  probability  is  set  to  1.0  the  decision  rule 
can  be  considered  to  be  a  preferential  policy,  i.e.,  tl  can  fire  only  when  it  is  not  in 
conflict  with  t2  (equivalent  to  saying  that  tl  can  fire  only  when  p3  has  no  token). 
It  is  clear  from  this  example  that  the  dynamics  of  some  Petri  nets  cannot  be  made 
explicit  without  decision  rules.  Once  the  static  properties  of  a  Petri  net,  i.e.,  the  net 
topology  and  initial  token  placement,  are  defined,  transitions  that  may  potentially 
be  in  conflict  with  each  other  can  be  identified.  A  set  of  such  transitions  is  called  a 
conflict  set.  Given  a  Petri  net  topology  and  initial  token  placement,  it  is  important 
to  identify  all  conflict  sets  and  assign  a  decision  rule  for  each  set.  See  [7]  and  [8]  for 
more  detailed  discussion  on  net  topology  and  conflict  sets. 

To  summarize,  the  dynamics  of  a  Petri  net  depends  on  specification  of  the  net 
structure  (topology),  the  initial  placement  of  the  tokens,  and  the  decision  rules  (if 
there  are  conflict  sets). 

3.2.  Stochastic  Petri  Nets 

The  original  Petri  nets,  as  described  above,  are  defined  without  the  notion  of  time; 
however,  in  timed  Petri  nets  some  explicitly  defined  timing  parameters  (“process- 
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ing  times”)  influence  the  evolution  of  the  state  of  the  nets.  Stochastic  Petn  nets 
(SPN’s)  are  those  timed  Petri  nets  in  which  the  processing  times  are  specified  prob¬ 
abilistically  via  distribution  functions.  Indeed,  the  relation  between  SPN’s  and  the 
original,  untimed  Petri  nets  is  analogous  to  that  between  semi-Markov  chains  and 
Markov  chains.  Different  versions  of  SPN’s  exist.  For  example,  SPN’s  can  be  di¬ 
vided  into  two  groups  depending  on  whether  the  processing  times  are  associated 
with  transitions  or  places8.  This  issue  is  particularly  important  in  modeling  of  car¬ 
diac  arrhythmias,  and  we  discuss  it  in  detail  in  this  subsection.  Another  source  of 
diversification  of  SPN’s  is  the  relationship  between  clocking  of  the  processing  times 
and  firing  of  the  tokens.  Reference  [6]  describes  various  “firing  policies”  that  result 
because  of  different  interpretations  of  this  relationship.  The  firing  policy  of  interest 
in  this  paper  is  called  the  “race”  policy  which  is  assumed  in  the  descriptions  below. 
Tfransition-timed  SPN.  In  transition-timed  SPN’s,  the  processing  times  are  as¬ 
sociated  with  transitions.  When  a  transition  is  enabled,  a  sample  of  the  random 
variable  representing  the  corresponding  processing  time  is  chosen,  and  the  transi¬ 
tion  must  wait  for  this  amount  of  time  before  it  can  fire.  It  is  possible  that  the 
transition  becomes  disabled  while  waiting  to  fire.  For  example,  in  Figure  10,  tO 
can  be  disabled  while  waiting  (regardless  of  decision  rule)  if  tl  becomes  enabled  and 
fires.  Note  that  in  the  figure  a  double  bar  represents  a  transition  whose  processing 
time  is  non-zero.  A  single  bar  represents  a  transition  with  zero  processing  time;  it 
fires  as  soon  as  it  is  enabled. 

Place-timed  SPN.  In  place-timed  SPN’s,  the  processing  times  are  associated  with 
places.  When  a  token  enters  a  place,  it  initially  becomes  “unavailable”  to  the  rest  of 
the  system  until  the  corresponding  probabilistically  chosen  processing  time  elapses. 
A  transition  is  not  enabled  unless  all  its  input  arcs  are  connected  to  places  with 
“available”  tokens.  Figure  11  shows  an  example  of  place-timed  SPN’s.  A  double 
circle  in  the  figure  represents  a  place  with  non-zero  processing  time,  while  a  single 
circle  represents  a  place  with  zero  processing  time. 

A  place-timed  SPN  can  be  converted  to  a  transition-timed  SPN  that  has  the 
same  dynamic  property;  however,  the  converse  of  this  statement  is  not  true  in  gen- 
*No  SPN  developed  so  far  has  attempted  to  assign  transition  times  to  both  transitions  and  places. 
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eral.  For  example,  the  transition-timed  SPN  fragment  depicted  in  Figure  10  has 
no  equivelence  in  the  place-timed  format9.  Nevertheless,  in  many  modeling  prob¬ 
lems  place-timed  SPN’s  offer  more  concise  structure  than  transition-timed  SPN’s 
(Fig.  12),  making  model  interpretation  easier.  This  is  especially  true  in  modeling 
of  the  cardiac  electrical  conduction  system,  and  the  issue  of  choosing  one  of  the  two 
SPN  formats  is  addressed  breifly  in  the  next  subsection.  Both  formats  can  be  seen 
regularly  in  the  literature,  although  transition-timed  SPN’s  seem  to  be  a  little  more 
popular. 

3.3.  SPN  for  Modeling  of  Cardiac  Conduction  System 

SPN  modeling  of  the  cardiac  electrical  conduction  system  poses  a  tradeoff  between 
using  the  place-timed  and  transition-timed  versions:  (VjPlace-timed  SPN’s  offer  a 
graphically  more  concise  format  to  represent  the  cardiac  system  than  transition- 
timed  SPN’s.  For  example,  compare  the  normal  heart  models  using  the  two  SPN 
formats(Fig.  2c  and  2d).  (Yi/Transition-timed  SPN’s  can  model  a  wider  range  of 
dynamics  than  place-timed  SPN’s.  Unfortunately,  some  dynamic  aspects  of  the  car¬ 
diac  system  can  only  be  modeled  with  the  transition-timed  format.  An  example  of 
such  is  activation  of  a  tissue  during  its  relative  refractory  period.  The  SPN  model 
of  relative  refractory  period  (Fig.  14d-l)  is  almost  identical  to  the  transition-timed 
SPN  of  Figure  10,  which  cannot  be  translated  into  a  place-timed  SPN.  We  em¬ 
phasize  that  the  graphical  conciseness  offered  by  place-timed  SPN’s  is  important, 
as  we  will  see  in  the  next  section  that  some  relatively  simple  models  of  cardiac 
arrhythmia  can  have  complex  topologies.  Graphical  conciseness  facilitates  not  only 
interpretation  of  the  models  but  also  representation  of  more  complex  physiological 
mechanisms.  On  the  other  hand,  the  transition-timed  format  is  necessary  to  repre¬ 
sent  certain  essential  dynamical  properties  of  the  cardiac  system.  This  dilemma  can 
be  avoided  in  the  following  way:  As  we  have  mentioned,  every  place-timed  SPN  has 
an  equivalent  transition-timed  SPN.  We  can  thus  graphically  represent  a  system 

®This  statement  holds  if  the  "race  firing  policy"  [6]  is  adopted.  In  the  literature  this  policy 
is  seldom  seen,  and  statements  asserting  equivalence  between  place-  and  transition-timed  formats 
can  be  found.  There  is,  however,  no  single  SPN  format  preferred  over  the  others  today,  and  such 
statements  are  not  well  founded. 
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with  a  place-timed  SPN  but  deal  with  its  dynamics  in  terms  of  its  transition-timed 
equivalence.  This  relation  between  place-timed  SPN’s  and  transition-timed  SPN’s 
is  analogous  to  that  between  high  level  computer  languages  and  assembly  codes  as 
compilers  translate  the  high  level  languages  into  assembly  codes  before  executable 
codes  are  generated. 

To  represent  a  model  with  a  place-timed  SPN  while  dealing  with  its  dynamics 
with  the  corresponding  transition-timed  SPN,  we  need  to  allow  the  place-time  for¬ 
mat  to  represent  just  as  wide  range  of  dynamics  as  the  transition-timed  format. 
We,  thus,  make  a  modification  in  the  place-timed  SPN  format  by  introducing  “in- 
terruptable  processing  times.”  A  place  with  an  interruptable  processing  time  has 
a  special  output  arc  through  which  a  token  inside  it  is  always  “available”  to  the  rest 
of  the  system.  Figure  13a  shows  such  a  place.  It  has  a  special  output  arc,  marked 
by  a  small  circle,  through  which  the  token  inside  it  is  always  available  to  the  rest 
of  the  system.  Thus,  while  the  transition  tO  cannot  fire  until  the  processing  time  at 
the  place  elapses,  tl  can  fire  as  soon  as  it  is  enabled.  Firing  of  tl,  therefore,  may 
“interrupt”  the  processing  time  assigned  to  the  place.  As  previously  mentioned, 
a  standard  place-timed  format  cannot  replace  the  transition-timed  SPN  fragment 
of  Figure  10,  but  with  an  interruptable  processing  time  that  is  possible(Fig.  13c). 
This  illustrates  the  usefulness  of  the  new  place-timed  SPN’s.  Although  this  modified 
place-timed  format  (i.e.,  using  interruptable  processing  times)  is  not  found  in  the 
literature,  its  execution  is  exactly  the  same  as  the  standard  transition-timed  format 
under  the  “race”  firing  policy.  For  example,  the  modified  place-timed  SPN  of  Fig¬ 
ure  13a  behaves  exactly  the  same  as  the  transition-timed  SPN  of  Figure  13b.  That 
is,  a  modified  place-timed  SPN  can  always  be  converted  into  a  standard  transition- 
timed  SPN,  whose  execution  policy  is  well-defined.  We  will  be  using  this  modified 
place-timed  SPN  format  in  the  next  section. 
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4.  SPN  Models  of  Cardiac  Arrhythmias 

In  Section  2  we  modeled  several  common  cardiac  arrhythmias  with  networks  of 
rhythm  and  transmission  elements.  In  this  section  we  implement  various  rhythm 
and  transmission  elements  with  SPN’s.  We  then  discuss  briefly  how  to  connect 
these  SPN  elements  to  form  SPN  cardiac  arrhythmia  models.  Finally,  we  present 
several  examples  of  such  SPN  models  of  cardiac  arrhythmias  at  the  end  of  the 
section.  The  graphical  notations  in  the  figures  of  this  section  are  the  same  as  the 
previous  sections.  In  particular,  the  modified  place-timed  SPN’s  are  denoted  as 
in  Figure  13:  double  circles  represent  places  with  non-zero  processing  times,  and 
output  arcs  marked  by  small  circles  are  the  “interruptable”  outlets  of  tokens  from 
such  places. 

4.1.  Fundamental  Building  Blocks 

Recall  form  Section  2  that  variations  in  autorhythmic  rates,  conduction  delays,  and 
refractory  periods  characterize  functions  of  individual  rhythm  and  transmission  ele¬ 
ments  which  are  in  turn  responsible  for  most  cardiac  arrhythmias.  We  now  present 
SPN  implementations  of  sub-elements  (“units")  controlling  these  timing  quantities. 
1) Autorhythmic  unit. 

Figure  14a  shows  the  basic  SPN  building  block  for  an  autorhythmic  unit.  The  pe¬ 
riod  of  the  rhythm  is  represented  by  the  processing  time  associated  with  the  single 
place  pO.  The  transition  tl  represents  the  output  of  the  autorhythmic  unit.  The 
incomplete  arcs  to  and  from  tl  are  parts  of  the  SPN  fragment  representing  the 
neighboring  tissue  which  receives  excitation  from  the  autorhythmic  unit.  The  re¬ 
ception  of  excitation  is  accomplished  by  firing  of  tl.  As  soon  as  the  token  becomes 
availabe  in  pO,  it  is  fired  back  into  pO  via  either  tO  or  tl  to  recycle  the  process. 
An  autorhythmic  activity  may  or  may  not  induce  activity  in  the  neighboring  tissue 
depending  on  whether  the  tissue  is  in  refractory  period  or  not.  But  if  the  tissue 
is  ready  to  be  excited,  the  activity  is  always  transferred  to  it.  Thus,  the  decision 
rule  for  the  conflict  set  {tO,  tl}  is  assigned  so  that  it  chooses  tl  preferentially  over  tO. 
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2) Conduction  units. 

A  conduction  unit  receives  excitation  from  one  tissue  (the  input  tissue),  waits  for  a 
probabilistically  specified  amount  of  time  (the  conduction  delay),  and  transfers  the 
excitation  to  another  tissue  (the  output  tissue).  Figure  14b  shows  the  SPN  imple¬ 
mentation  of  the  basic  conduction  unit.  If  the  token  in  the  place  pi  is  “available” 
(meaning  that  the  conduction  tissue  has  come  out  of  the  refractory  period  and  is 
excitable),  the  transition  tO  fires  when  it  receives  excitation  from  the  input  tissue. 
The  token  then  enters  the  place  pO  whose  processing  time  represents  the  conduction 
delay.  After  the  delay,  the  token  is  fired  back  into  pi  by  tl  or  t2.  When  the  token  is 
in  pO,  tl  is  always  enabled  while  t2  is  enabled  only  when  the  output  tissue  is  ready 
to  receive  the  excitation.  Since  excitation  is  always  transferred  to  output  tissue 
not  in  its  refractory  period,  the  decision  rule  is  assigned  so  that  t2  is  preferentially 
chosen  over  tl  whenever  they  are  in  conflict.  The  processing  time  associated  with 
pi  represents  the  time  required  for  the  conduction  unit  to  become  ready  to  receive 
a  new  excitation  after  it  has  processed  a  preceding  excitation;  thus,  this  processing 
time  can  be  considered  as  a  collective  absolute  refractory  period  for  all  the  cells 
in  the  conduction  tissue.  The  roles  of  relative  refractory  periods  in  conduction  of 
excitations  is  discussed  later. 

3)  Refractory  periods. 

There  are  two  kinds  of  refractory  periods-absolute  and  relative.  When  the  tissue  is 
in  an  absolute  refractory  period,  an  oncoming  excitation  is  blocked.  On  the  other 
hand,  when  the  tissue  is  in  a  relative  refractory  period,  an  oncoming  excitation 
induces  the  tissue  to  be  activated  with  abnormal  characteristics  such  as  elongated 
conduction  delay  (e.g.,  Wenckebach).  Figure  14c  shows  a  basic  absolute  refractory 
unit.  The  processing  time  assigned  to  the  single  place  is  the  absolute  refractory  pe¬ 
riod,  and  the  transition  cannot  fire  unless  the  token  in  the  place  becomes  available. 
Figure  14d-2  shows  a  basic  relative  refractory  unit,  using  a  place  with  an  interrupt- 
able  processing  time  introduced  in  Section  3.  The  processing  time  associated  with 
the  place  represents  the  relative  refractory  period.  The  firing  of  tO  represents  the 
normal  course  of  action  where  the  tissue  becomes  completely  ready  to  receive  a  new 
excitation.  On  the  other  hand,  the  firing  of  tl  represents  a  premature  activation 
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of  the  tissue  by  an  external  excitation  (an  interrupt),  and  the  tissue  is  expected  to 
display  some  abnormal  activities.  Figure  14d-l  shows  the  transition-timed  equiva¬ 
lent  of  the  relative  refractory  unit.  As  mentioned  in  the  previous  section,  these  two 
SPN’s  display  exactly  the  same  dynamics. 

4.2.  SPN  Model  of  the  Normal  Heart 

Using  the  building  blocks  developed  above,  a  SPN  model  of  the  normal  heart  can 
be  implemented.  We  have  shown  in  Section  2  that  the  normal  heart  can  be  modeled 
with  three  elements  —  the  SA  node,  AV  node,  and  ventricles  (Fig.  2b),  and  these 
elements  can  respectively  implemented  as  the  autorhythmic  unit,  conduction  unit, 
and  absolute  refractory  unit  described  above.  The  resulting  SPN  model  is  shown  in 
Figure  2c.  The  output  of  the  model  is  the  timing  of  generation  of  P  and  R  waves,  and 
these  are  represented  by  firing  of  the  transitions  tO,  tl,  and  t3  —  tO  and  tl  for  P  wave 
and  t3  for  R  wave.  The  parameters  of  the  model  are  the  processing  times  associated 
with  the  places  pO,  pi,  p2,  and  p3,  and  these  correspond  with  the  autorhythmic 
interval  of  the  SA  node,  the  refractory  period  of  the  AV  node,  conduction  delay 
of  the  AV  node,  and  the  refractory  period  of  the  ventricles,  respectively.  The 
transition-timed  equivalent  is  shown  in  Figure  2d,  and  note  that  the  model  in  this 
format  requires  twice  as  many  nodes  as  the  model  in  the  place-timed  format.  Figure 
2e  shows  the  output  of  a  simulation  run  of  the  SPN  model  in  Figure  2c.  (As 
mentioned  before,  short  and  long  vertical  lines  represent  P  waves  and  R  waves, 
respectively,  and  tick  marks  have  one  second  intervals.)  The  processing  times  used 
in  the  simulation  run  are  listed  on  Table  1. 

The  simplicity  of  the  models  in  Figure  2  has  come  about  because  unnecessary 
physiological  details  are  deliberately  left  out.  For  example,  in  the  signal  flow  block 
diagram  model  (Fig.  2b),  uni-directional  transmission  element  rather  than  more 
complex  bi-directional  transmission  element  was  used  to  represent  the  AV  node 
because  it  is  known  that  no  retrograde  conduction  is  observed  in  a  normal  P,R  wave 
sequence.  The  rhythm  element  representing  the  SA  node  does  not  need  an  input 
terminal  also  because  of  the  absence  of  retrograde  conduction.  We  emphasize  that 
the  normal  heart  is  still  capable  of  conducting  excitations  in  the  retrograde  direction, 


but  modeling  at  such  level  of  physiological  accuracy  is  not  necessary  to  capture  the 
behavior  of  a  normal  P,R  wave  sequence  correctly.  Using  such  simpler  variations 
of  the  rhythm  and  transmission  elements  allows  simpler  SPN  implementations  of 
them.  In  this  case  we  were  able  to  directly  substitute  the  autorhythmic  unit  from 
Subsection  4.1  for  the  rhythm  element  representing  the  SA  node,  the  conduction 
unit  for  the  transmission  element  of  the  AV  node,  and  the  absolute  refractory  unit 
for  the  rhythm  element  of  the  ventricles.  We  have  seen  in  Section  2,  however, 
that  more  complex  variations  of  these  elements  are  needed  to  model  most  cardiac 
arrhythmias.  In  the  next  subsection  we  extend  the  basic  SPN  units  of  the  previous 
subsection  to  implement  a  wider  variety  of  rhythm  and  transmission  elements. 

4.3.  Rhythm  and  Transmission  Elements 

The  rhythm  and  transmission  elements  are  aggregations  of  the  fundamental  building 
blocks  developed  above.  SPN  implementations  of  various  rhythm  and  transmission 
elements  are  presented. 

Rhythm  Elements. 

The  most  simple  rhythm  elements  are  the  pacemaker  (Fig.  15a)  and  terminal 
(Fig.  15f).  An  example  of  a  pacemaker  is  the  SA  node  in  the  normal  heart;  the 
sole  function  of  the  element  is  to  generate  excitation  periodically.  The  single  pa¬ 
rameter  for  the  pacemaker  is  the  period  of  the  autorhythmic  excitation,  and  this  is 
represented  by  the  processing  time  at  the  place.  An  example  of  a  terminal  is  the 
ventricle  in  the  normal  heart;  the  element  simply  receives  and  absorbs  incoming  ex¬ 
citations.  The  single  parameter  of  the  terminal  is  the  absolute  refractory  period.  If 
the  terminal  depicted  on  Figure  15f  represents  the  ventricle,  firing  of  the  transition 
represents  the  generation  of  the  R  wave. 

In  reality,  pacemaker  tissues  can  receive  excitation  (for  example,  the  SA  node 
can  receive  retrograde-conducting  excitation)  and  can  be  reset  or  stunned.  A  more 
accurate  model  of  pacemakers  should  have  such  dynamic  properties.  Moreover,  a 
way  to  make  the  model  more  accurate  than  this  is  to  include  an  absolute  refractory 
period  during  which  the  pacemaker  tissue  cannot  be  reset  or  stunned.  These  func¬ 
tional  variations  of  the  rhythm  element  can  be  used  in  arrhythmia  modeling,  and 
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corresponding  SPN  building  blocks  are  presented  in  Figure  15: 
the  basic  pacemaker  (Fig.  15a) 
resettable  pacemaker  (Fig.  15b) 
stunnable  pacemaker  (Fig.  15c) 

resettable  pacemaker  with  refractory  period  (Fig.  15d) 
stunnable  pacemaker  with  refractory  period  (Fig.  15e) 
the  basic  terminal  (Fig.  15f) 

Transmission  Elements. 

The  most  simple  transmission  element  is  the  basic  uni-directional  transmission  ele¬ 
ment  depicted  in  Figure  16a  and  employed  in  the  construction  of  the  normal  heart 
model  (Fig.  2b).  One  of  the  functional  variations  is  the  transmission  element  with  a 
probability  of  conduction,  shown  in  Figure  16b,  which  is  necessary  in  the  modeling 
of  second  degree  AV  block.  An  extra  transition  tO  has  been  added.  The  transitions 
tO  and  tl  are  in  conflict;  firing  of  tO  does  not  activate  the  transmission  element  while 
firing  of  tl  does.  The  decision  rule  for  the  conflict  set  {tO,  tl}  is  a  probabilistic 
policy  where  the  probability  “q”  of  choosing  tl  is  equal  to  the  probability  of  the 
conduction. 

Another  functional  variation  is  the  bi-directional  transmission  element.  It  is 
more  realistic  than  the  uni-directional  element,  while  its  implementation  is  slightly 
more  complex  than  a  simple  combination  of  two  uni-directional  elements.  The  com¬ 
plexity  is  due  to  the  necessity  of  modeling  annihilation  of  two  excitations  colliding 
inside  the  element.  Figure  16c  shows  the  bi-directional  transmission  element.  It 
is  basically  two  uni-directional  elements  superimposed,  with  the  addition  of  the 
transition  t5  and  the  use  of  interruptable  processing  times  at  pO  and  pi.  When 
excitations  collide  in  the  element,  t5  fires,  and  no  activity  reaches  the  outputs. 

Another  functional  variation  is  the  transmission  elements  capable  of  generating 
the  Wenckebach  phenomenon.  As  discussed  in  Section  2  the  Wenckebach  phe¬ 
nomenon  is  charcterized  by  an  incrementally  increasing  P,R  sequence.  An  abnor¬ 
mally  long  refractory  period  at  the  AV  node  is  thought  to  be  responsible  for  the 
phenomenon.  Specifically,  after  the  AV  node  is  excited  by  a  normal  excitation  from 
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the  SA  node,  the  relative  refractory  period  of  the  node  is  just  long  enough  such 
that  the  next  normal  excitation  arrives  at  the  AV  node  during  a  late  part  of  the 
relative  refractory  period.  The  AV  node  still  becomes  excited,  but  the  subsequent 
conduction  delay  as  well  as  absolute  and  relative  refractory  periods  are  longer  than 
previous  ones.  Consequently,  the  next  normal  excitation  tends  to  arrive  at  the  AV 
node  during  an  even  earlier  part  of  the  relative  refractory  period,  resulting  in  yet 
longer  conduction  delays  and  refractory  periods.  Such  an  incremental  increase  in 
the  AV  node  conduction  delay  (and  thus  the  P,R  wave  interval)  continues  until 
the  normal  excitation  arrives  at  the  node  during  an  absolute  refractory  period,  re¬ 
sulting  in  a  blocked  conduction  (and  thus  a  missing  R  wave).  When  a  conduction 
block  occurs,  the  AV  node  gains  enough  time  to  recover  from  its  refractory  period 
completely;  then,  the  process  described  above  begins  again. 

Figure  16d  shows  an  SPN  implementation  of  a  transmission  element  that  cap¬ 
tures  such  cyclic  process.  The  relative  refractory  period  is  represented  by  the  inter- 
ruptable  processing  time  of  pO.  The  normal  course  of  token  movement  is  pO  -  pi  - 
p2  -  p3  -  pO.  However,  when  the  element  is  excited  during  its  relative  refractory 
period,  the  token  travels  the  abnormal  course  pO  -  p4  -  p5  -  pO.  The  conduction 
delays  through  the  element  are  represented  by  p2  and  p4,  while  the  absolute  refrac¬ 
tory  periods  are  represented  by  p3  and  p5.  For  the  Wenckebach  phenomenon,  the 
processing  times  for  p4  should  be  set  longer  than  for  p2  while  that  for  p3  should  be 
longer  than  for  p5.  We  have  mentioned  in  Section  2  that  this  element  should  have 
“memory.”  The  element  in  this  example  registers  whether  the  previous  excitation 
has  arrived  during  its  relative  refractory  period  or  not  by  the  location  of  the  token, 
i.e.,  whether  the  token  is  in  the  abnormal  loop  or  the  normal  loop.  Note  that  the 
token  takes  a  longer  time  to  go  through  the  abnormal  loop  so  that  the  chance  of 
receiving  the  next  external  excitation  while  it  is  in  absolute  or  relative  refractory 
period  is  greater.  To  model  common  Wenckebach  conditions,  more  than  one  abnor¬ 
mal  loop  may  be  necessary.  Such  a  Wenckebach  model  is  presented  in  Subsection 
4.5. 

Several  functional  variations  of  the  transmission  element  have  been  discussed 
here.  Although  their  coverage  is  not  exhaustive,  most  of  common  cardiac  arrhyth- 
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mias  can  be  modeled  using  the  variations  of  the  elements  discussed  above.  The 
SPN  implementations  of  these  are  presented  in  Figure  16: 

the  basic  uni-directional  transmission  element  (Fig.  16a) 
uni-directional  transmission  element  with  probability  of  block  (Fig.  16b) 
bi-directional  transmission  element  (Fig.  16c) 
a  stunnable  transmission  element  (Fig.  16d) 

4.4.  Interfacing  the  Elements 

Here,  we  discuss  ways  to  connect  the  SPN  representations  of  rhythm  and  transmis¬ 
sion  elements  described  above. 

Case  I:  connecting  a  single  output  to  multiple  inputs.  Figure  17a  shows  a  rhythm 
element  sending  excitation  to  two  transmission  elements.  When  the  rhythm  element 
sends  an  excitation,  each  of  the  two  receivers  (the  transmission  elements)  can  be 
in  two  states —  ready  to  receive  the  excitation  or  not  ready  to  receive.  Thus,  the 
number  of  states  of  the  receivers  is  four  (=  22),  and  four  transitions  are  used  to 
implement  this  interface  with  SPN  (Fig.  17b).  Each  of  the  four  transitions  represent 
the  following  cases: 

tO  -  neither  of  the  transmission  element  is  ready  to  receive, 
tl  -  element  A  is  ready  to  receive,  but  not  element  B. 
t2  -  both  elements  are  ready  to  receive. 
t3  -  element  B  is  ready  to  receive,  but  not  element  A. 

In  general,  when  a  single  output  is  trying  to  distribute  an  excitation  to  N  inputs, 
2n  transitions  are  required. 

Case  II:  connecting  multiple  outputs  to  a  single  input.  Figure  18a  shows  three 
rhythm  elements  sending  excitation  to  a  transmission  element.  Since  each  sender 
may  or  may  not  be  able  to  excite  the  receiver,  six  (=  2  x  3)  transitions  are  needed 
(Fig.  18b).  In  general,  when  N  outputs  are  trying  to  access  to  a  single  input,  2 N 
transitions  are  required. 
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4.5.  SPN  Models  of  Cardiac  Arrhythmias 

SPN  implementations  of  the  cardiac  arrhythmias  discussed  in  Section  2  are  pre¬ 
sented  in  Figures  3c  to  8c.  The  results  of  simulation  based  on  these  SPN’s  are 
presented  in  Figures  3b  to  8b,  in  which  the  normal  P  waves  are  represented  by  short 
lines  and  normal  R  waves  by  long  lines.  The  values  of  the  model  parameters  used 
in  the  simulation  runs  are  listed  in  Table  1. 

Atrial  Premature  Beat  (Fig.  3c) 

The  block  diagram  model  of  Figure  3a  is  implemented  with  an  SPN.  Two  “reset¬ 
table  pacemakers”  are  used  to  represent  the  normal  (SA  nodal)  and  ectopic  rhythm 
sources  in  the  atria.  Firing  of  one  source  resets  the  other  source;  thus,  after  every 
beat  the  two  sources  race  to  initiate  the  next  excitation.  The  processing  times  at 
pO  and  pi  correspond  with  the  autorhythmic  intervals  of  the  ectopic  and  normal 
sources,  respectively.  The  probability  density  function  of  the  processing  time  at  pO 
has  a  broader  distribution  than  the  one  at  pi,  so  that  with  a  certain  probability 
the  ectopic  source  can  generate  excitation  at  a  noticablly  shorter  interval  than  the 
normal  source.  The  firing  of  tO  or  tl  represents  an  excitation  of  the  ectopic  source, 
and  it  produces  a  P  wave  with  an  abnormal  morphology,  denoted  by  a  “P”  with  an 
overbar  in  the  figure.  The  firing  of  t2  or  t3  represents  an  SA  nodal  excitation  pro¬ 
ducing  a  normal  P  wave.  The  rest  of  the  model  is  similar  to  the  normal  heart  model. 
The  places  p3  and  p2  represent  the  AV  nodal  conduction  delay  and  its  refractory 
period,  respectively,  and  they  are  parts  of  the  “basic  uni-directional  transmission 
element.”  The  ventricles  are  modeled  by  the  “terminal”  element  (p 4  and  t5).  Firing 
of  t5  produces  a  R  wave.  In  the  simulation  result  (Fig.  3b),  the  short  vertical  lines 
with  small  squares  at  their  tips  represent  P  waves  initiated  by  the  ectopic  source 
(i.e.,  the  premature  beats). 

Ventricular  Premature  Beat  (Fig.  4c) 

This  model,  whose  block  diagram  was  introduced  in  Fig.  4a,  consists  of  two  “reset¬ 
table  pacemakers  with  refractory  periods”  and  a  single  “bi-directional  transmission 
element.”  The  places  pi  and  pO  are  the  autorhythmic  interval  and  absolute  refrac¬ 
tory  period  at  the  SA  node,  respectively,  while  p7  and  p6,  respectively,  are  those 
at  the  ventricles.  The  conduction  delay  and  absolute  refractory  period  at  the  AV 
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node  are  represented  by  p3  and  p4,  respectively,  for  the  antegrade  (norma!)  direc¬ 
tion  and  by  p2  and  p5,  respectively,  for  the  retrograde  direction.  Correspondence 
between  firing  of  the  transitions  and  production  of  the  waves  are  as  follows:  t2 
and  t3  for  normal  P,  t5  for  normal  R,  t6  and  t7  for  ectopic  (premature)  R,  and  tl 
for  retrograde  P.  In  the  simulation  result(Fig.  4b),  the  R  waves  and  P  waves  with 
small  square  on  the  tips  denote  the  ectopic  R  waves  and  retrograde  P  waves,  respec¬ 
tively.  Observe  that  a  pairing  of  a  normal  P  wave  and  an  ectopic  R  wave  occuring 
at  about  the  same  time  is  followed  by  neither  a  normal  R  wave  nor  a  retrograde 
P  wave.  Nearly  simultaneous  autorhythmic  excitations  at  the  two  rhythm  sources 
of  the  model  produce  a  normal  P  and  ectopic  R  waves,  but  the  resulting  flows  of 
excitations  have  collided  inside  the  transmission  element  and  blocked  each  other. 
Bigeminy  (Fig.  5c) 

As  described  in  Section  2,  we  use  the  concept  of  “reentrant  pathway”  to  explain 
the  mechanism  of  bigeminy.  The  model  is  almost  the  same  as  the  normal  heart 
model.  The  only  difference  is  the  presence  of  the  reentrant  pathway  represented  by 
the  basic  uni-directional  transmission  element  consisting  of  p3  and  p5.  When  the 
excitation  is  passed  from  the  AV  node  to  the  ventricles  (which  is  modeled  by  the 
“terminal”  element  consisting  of  p4)  via  the  transition  t4,  the  reentrant  pathway 
is  excited.  The  reentrant  pathway  delays  the  excitation  by  an  amount  specified  by 
the  processing  time  at  p3  and  then  passes  it  to  the  ventricles  via  t5.  Firing  of  t5 
produces  abnormal  R  waves  which  are  denoted  by  long  lines  with  small  squares  in 
the  simulation  result(Fig.  5b).  The  time  intervals  between  a  normal  R  wave  and 
the  succeeding  abnormal  R  wave  are  fairly  constant,  exhibiting  a  coupling  interval. 
Note  that  a  normal  P  wave  follows  every  abnormal  R  wave.  In  reality  such  a  P 
wave  cannot  be  easily  discerned  in  an  ECG  tracing  because  of  the  typical,  broad 
shape  of  the  abnormal  R  waves  of  bigeminy. 

Second  Degree  AV  Block  (Fig.  6c) 

The  topology  of  this  model  is  exactly  the  same  as  that  of  the  normal  heart,  but  in 
this  model  “uni-directional  trasmission  element  with  probability  of  conduction”  is 
used  to  represent  the  AV  node.  Thus,  the  decision  rule  for  the  conflict  set  {t0,  tl}  is 
no  longer  determinstic  (i.e.,  to  preferentially  choose  tl  over  tO).  Whenever  a  conflict 
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exists,  tl  is  chosen  over  tO  with  a  probability  equal  to  the  probability  of  conduction. 
Figure  6b  is  the  result  of  simulation  of  a  “4:3  Block,”  where  the  probability  of 
conduction  is  Since  in  this  model  conduction  blocks  are  independent  events, 
they  can  occur  in  consecutive  beats.  In  reality,  however,  such  consecutive  blocks 
are  rare.  Making  sure  that  conduction  blocks  are  isolated  can  be  accomplished  at 
the  expense  of  a  more  complex  model  with  more  memory.  Specifically,  decision  by 
the  model  to  cause  a  block  in  a  particular  beat  must  be  influenced  by  whether  or 
not  a  block  has  occurred  in  the  previous  beat,  i.e.,  the  model  must  have  a  “memory” 
of  an  event  in  the  previous  beat.  Although  such  a  model  is  not  presented  in  this 
paper,  the  model  of  the  Wenckebach  phenomenon  described  next  contains  a  memory 
feature  of  the  same  general  type  as  that  would  be  needed  to  obtain  a  more  accurate 
model  of  second  degree  AV  block. 

Wenckebach  (Fig.  7c) 

The  block  diagram  of  this  Wenckebach  model  consists  of  two  rhythm  elements 
and  a  transmission  element(Fig.  7a).  Two  “basic  rhythm  elements”  are  used  to 
represent  the  SA  node  (pO)  and  ventricles  (pl2).  The  single  transmission  element  is 
a  “stunnable  transmission  element”  similar  to  the  one  described  in  Subsection  4.3 
and  Figure  16d,  but  this  element  has  two  extra  “abnormal  loops.”  The  loops,  in 
the  order  of  increasing  transit  times,  are:  (i)pl  -  p2  -  p3  -  p4  -  p8,  which  is  the 
normal  loop,  (ii)pl  -  p2  -  p3  -  p5  -  p9,  (iii)pl  -  p2  -  p6  -  plO,  and  (iv)pl  - 
p7  -  pll.  The  processing  times  at  p4,  p5,  p6,  and  p7  (in  the  order  of  increasing 
length)  represent  the  AV  conduction  delays  (thus  the  P,R  intervals).  Those  at  p8, 
p9,  plO,  and  pll  are  the  absolute  refractory  periods  of  the  AV  node.  The  relative 
refractory  period  is  divided  into  three  sub-periods  of  equal  lengths  —  early  (pi), 
middle  (p2),  and  late(p3).  Each  of  these  sub-periods  is  assigned  an  interruptable 
processing  time.  An  interrupt  occurs  when  the  SA  node  (pO)  excites  while  the  token 
of  the  AV  node  is  in  one  of  these  relative  refractory  sub-periods,  initiating  an  entry 
into  the  corresponding  abnormal  loop.  P  waves  are  produced  by  the  firing  of  to, 
tl,  t2,  t3,  or  t4,  while  R  waves  by  the  firing  of  t8,  tlO,  1 12,  or  tl4.  The  simulation 
result(Fig.  7b)  clearly  displays  the  Wenckebach  phenomenon. 


Wolfe-Parkinson- White  Syndrome  (Fig.  8c) 

This  model  uses  “basic”  elements  for  the  SA  node  (pO-autorhythmic  interval)  and 
AV  node  (p3-conduction  delay,  p2-refractory  period).  The  ventricles  are  divided 
into  two  parts:  one  part  receives  excitation  from  the  bundle  of  Kent  and  produces 
delta  waves,  while  the  other  part  receives  excitation  normally  and  produces  normal 
R  waves.  Both  of  these  are  modeled  with  the  “terminal”  rhythm  elements  (p5  and 
p6).  The  bundle  of  Kent  is  represented  by  the  uni-directional  transmission  element 
consisting  of  p4  (conduction  delay)  and  pi  (refractory  period).  Delta  waves  and 
normal  R  waves  are  produced  when  t6  and  t5  fire,  respectively.  P  waves  are  pro¬ 
duced  when  tO,  tl,  t2,  or  t3  fires.  In  the  simulation  result(Fig.  8b)  the  delta  waves 
are  represented  by  short  lines  with  squares  at  the  tips,  and  they  are  immediatedly 
followed  by  R  waves  as  they  should.  Note  that  in  the  model  the  two  parts  of  the 
ventricles  (p5  and  p6)  could  have  been  connected  with  a  bi-directional  transmission 
elements  since  physiologically  they  are  parts  of  a  single  block  of  tissue.  It  is  clear 
from  the  simulation  result,  however,  that  these  two  parts  are  excited  almost  simul¬ 
taneously  so  that  collision  of  excitations  in  such  a  transmission  element  is  evident. 
Thus,  to  achieve  the  most  concise  modeling  of  the  P,R  sequence  such  an  element  is 
not  included. 

This  section  along  with  Section  2  has  illustrated  a  systematic  way  in  which 
we  can  derive  SPN’s  for  cardiac  arrhythmia  models.  Specifically,  the  physiological 
mechanism  of  an  arrhythmia  is  first  described  as  a  signal  flow  block  diagram  with 
rhythm  and  transmission  elements.  The  block  diagram  may  be  simplified  greatly 
by  removing  physiological  details  unnecessary  to  characterize  the  P,R  sequence  of 
the  particular  arrhythmia  under  study.  Such  a  high  level  description  can  then  be 
translated  into  an  SPN  in  an  element-by-element  fashion  using  the  implementations 
presented  in  this  section  (and  straightforward  extensions  of  them  if  needed).  Thus, 
we  have  a  method  to  derive  the  SPN  structure  (topology),  initial  token  placement, 
and  decision  rules.  The  last  piece  of  information  required  to  complete  the  model 
are  the  parameters,  i.e.,  probability  distributions  for  the  processing  times  and  prob¬ 
abilistic  descision  rules.  We  emphasize  that  two  sets  of  parameters  assigned  to  the 
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same  SPN  structure  can  lead  to  two  outputs  drastically  different  from  each  other. 
For  example,  the  SPN  structure  for  Bigeminy  (Fig.  5c)  can  be  programmed  to  out¬ 
put  a  P,R  sequence  of  Trigeminy,  a  different  arrhythmia  class10,  just  by  increasing 
the  processing  time  at  the  place  p3  by  one  second. 

Given  a  P,R  sequence  to  be  modelled,  some  parameters  such  as  the  SA  nodal 
rhythms  and  AV  conduction  delays  have  direct  correspondence  to  the  given  data  and 
are  thus  relatively  easy  to  assign;  however,  many  other  parameters  such  as  refractory 
periods  cannot  be  directly  observed  from  the  data  and  must  be  chosen  intelligently. 
The  parameter  values  listed  in  Table  1  are  chosen  so  that  the  simulation  results 
match  with  the  qualitative  characteristics  of  a  P,R  wave  sequence  for  a  given  cardiac 
arrhythmia  as  described  in  Section  2,  under  the  constraint  that  each  value  chosen 
must  be  physiologically  plausible. 


10An  abnormal  R  wave  appears  every  three  cycle*  in  TVigeminy  in  contrast  to  every  two  cycles  in 
Bigeminy. 


5.  Conclusion 


In  this  paper  we  have  described  a  systematic  procedure  for  constructing  concise 
mathematical  models  capturing  the  cardiac  mechanisms  responsible  for  the  P,R 
wave  patterns  characteristic  of  different  arrhythmias.  In  particular,  there  is  a 
straightforward  way  in  which  one  proceeds  from  a  phenomenological  descriptions  to 
a  block  diagram  and  then  to  a  stochastic  Petri  net.  Since  the  methodology  employs 
characterizations  of  some  fundamental  properties  of  electrically  active  tissues  in  the 
heart,  it  is  more  robust  than  methods  taking  direct,  “black-box”  approaches  to 
describe  the  wave  sequences.  Several  important  areas  for  further  work  can  be  iden¬ 
tified.  First,  as  pointed  out  at  the  end  of  the  preceding  section,  there  is  the  issue  of 
how  to  select  the  parameter  values  systematically  to  realize  a  desired  P,R  sequence. 
The  development  of  an  automatic  algorithm  that  estimates  these  parameters  based 
on  a  given  P,R  sequence  and  a  set  of  physiological  constraints  (corresponding  to 
particular  cardiac  rhythms)  would  be  of  significant  value.  Such  an  algorithm  would 
also  be  needed  for  the  parameter  estimation  problem  in  a  cardiac  rhythm  classifi¬ 
cation  scheme,  as  each  model  must  be  “tuned”  to  the  observed  P,R  sequence  to  be 
classified.  Other  issues  related  with  rhythm  classification  are  wave  detection  and 
hypothesis  testing.  In  wave  detection,  one  can  imagine  integrating  the  physiological 
information  contained  in  the  model  with  the  process  of  identifying  wave  types  and 
timings  on  ECG  recordings.  In  hypothesis  testing,  an  efficient  method  to  evaluate 
“likelihoods”  for  alternate  hypothesized  SPN  models  must  be  developed. 
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Figure  1 

Building  blocks  for  the  signal  flow  models:  a  rhythm  element(a)  and  transmission 
element  (b). 


Figure  2 

Models  of  normal  heart  P,R  sequences:  (a)-a  physiologically  motivated  signal  flow 
model,  (b)-a  simplified  signal  flow  model,  (c)-SPN  implementation  of  (b)  using 
the  place-timed  format,  and  (d)-SPN  implementation  using  the  transition-timed 
format.  A  P,R  sequence  output  of  model  (c)  is  shown  in  (e). 

Figure  3 

Atrial  premature  beat:  signal  flow  model(a)  and  an  example  of  P,R  sequence(b) 
obtained  by  the  simulation  of  the  SPN  implementation (c)  of  the  model. 

Figure  4 

Ventricular  premature  beat:  signal  flow  model(a)  and  an  example  of  P,R  sequence(b) 
obtained  by  the  simulation  of  the  SPN  implementation (c)  of  the  model. 


Figure  5 

Bigeminy:  signal  flow  model(a)  and  an  example  of  P,R  sequence(b)  obtained  by  the 
simulation  of  the  SPN  implementation (c)  of  the  model. 


Figure  6 

Second  degree  AV  block:  signal  flow  model(a)  and  an  example  of  P,R  sequence(b) 
obtained  by  the  simulation  of  the  SPN  implementation(c)  of  the  model. 


Figure  7 

Wenckebach  phenomenon:  signal  flow  model  (a)  and  an  example  of  P,R  sequence(b) 
obtained  by  the  simulation  of  the  SPN  implementation (c)  of  the  model. 


Figure  8 

Wolfe-Parkinson- White  syndrome:  signal  flow  model(a)  and  an  example  of  P,R 
sequence(b)  obtained  by  the  simulation  of  the  SPN  implementation (c)  of  the  model. 

Figure  9 

An  example  of  Petri  net. 

Figure  10 

A  portion  of  a  transition-timed  SPN. 

Figure  11 

An  example  of  place-timed  SPN. 

Figure  12 

Equivalence  of  two  SPN  formats:  the  place-timed  SPN  fragment  in  (a)  can  be 
translated  into  the  transition-timed  format  as  shown  in  (b),  using  additional  nodes. 

Figure  13 

A  modified  place-timed  format  that  includes  interruptable  processing  time:  the  SPN 
fragment  in  the  modified  place-timed  format  shown  in  (a)  is  dynamically  equivalent 
with  that  in  the  transition-timed  format  shown  in  (b).  Also,  the  SPN  fragment  in  (c) 
(modified  place-timed)  is  dynamically  equivalent  with  that  in  Figure  10  (transition- 
timed),  which  cannot  be  translated  into  the  standard  place-timed  format  under  the 
“race”  execution  policy  defined  in  the  reference  [6]. 

Figure  14 

Fundamental  units  for  SPN  implementation  of  signal  flow  models. 

Figure  15 

SPN  implementations  of  common  rhythm  elements.  The  processing  time  assigned  to 
a  places  is  period  of  autorhythmicity  (denoted  as  “rhythm”  in  the  figure),  absolute 
refractory  period  (“refra”),  or  duration  of  stunning  (“stun”). 


Figure  16 

SPN  implementation  of  common  transmission  elements:  (a)-basic  uni-directional, 
(b)-uni-directional  with  probability  of  conduction  “q”,  (c)-bi-directional,  (d)-uni- 
directional  and  stunnable. 

Figure  17 

An  example  of  interfacing  SPN  blocks  —  fan-out:  (a)-The  rhythm  element  sends 
its  output  to  two  transmission  elements.  (b)-The  SPN  implementation  of  interface 
of  the  three  elements  in  (a). 

Figure  18 

An  example  of  interfacing  SPN  blocks  —  fan-in:  (a)-Each  of  the  three  rhythm 
elements  can  activate  the  transmission  element.  (b)-The  SPN  implementation  of 
interface  of  the  four  elements  in  (a). 


Second  Degree  AV  Block 


AV  node  (with 
Wenckebach  dynamics) 


Wenckebach 


(a)  Pacemaker 


(b)  Resettable  Pacemaker 


(d)  Resettable  Pacemaker 
with  absolute  refractory 
period 


(e)  Stunnable  Pacemaker 
with  absolute  refractory  period 


Table  1:  Model  Parameters 


Normal  Heart  (Fig.  2c) 
pO  (SA  rhythm) 
pi  (AV  refractory  period) 
p2  (AV  conduction  delay) 
p3  (Vent,  refractory  period) 

Atrial  Premature  Beat  (Fig.  3c) 
pO  (ectopic  rhythm) 

*non-uniform  distribution:  Prob(0.6) 
pi  (SA  rhythm) 
p2  (AV  refractory  period) 
p3  (AV  conduction  delay) 
p4  (Vent,  refractory  period) 

Ventricular  Premature  Beat  (Rg.  4c) 
pO  (SA  refractory  period) 
pi  (SA  rhythm) 

p2  (retrograde  conduction  delay) 
p3  (AV  conduction  delay) 
p4  (AV  refractory  -  normal) 
p5  (AV  refractory  -  retrograde) 
p6  (Ventricles  refractory) 
p7  (ectopic  rhythm) 

•non-uniform  distribution:  Prob(0.6) 

Bigeminy  (Fig.  5c) 
pO  (SA  rhythm) 
pi  (AV  refractory  period) 
p2  (AV  conduction  delay) 
p3  (reentrant  conduction  delay) 
p4  (Vent,  refractory  period) 
p5  (reentrant  refractory  period) 

Second  Dagisa  AY.  Bloch  (Fig.  6c) 
pO  (SA  rhythm) 
pi  (AV  refractory  period) 
p2  (AV  conduction  delay) 
p3  (Vent,  refractory  period) 

Decision  Rule  for  ( to,  tl } 


(0.92,  0.96.1.00,1.04.1.08} 
{0.3} 

{0.08,  0.10} 

{0-3} 


{0.6,  0.8, 1.5  }* 

Prob(0.8)  -  0.15;  Prob(1.5)  -  0.70. 

{0.92,  0.96,1.00,1.04,1.08} 
{0-3} 

{0.08,  0.10} 

{0.3} 


{0-3} 

{ 0.92,  0.96, 1.00, 1.04, 1.08 } 
{0.14,  0.16} 

{0.08,  0.10} 

{0.3} 

{0.3} 

{0.3} 

{0.6,  0.8, 1.5}* 

Prob(0.8)  -  0.15;  Prob(1.5)  -  0.70. 


{0.92,  0.96,1.00,1.04,1.08} 
{0.3} 

{0.08,  0.10} 

{  0.7,  0.8 } 

{0.5} 

{0} 


{0.92,  0.96,1.00,1.04,1.08} 
{0.3} 

{0.08,  0.10} 

{0.3} 

Prob(tO)  -  0.25;  Prob(t1)«0.75. 


^T.TWC'TT 


>:» 


# 

SP 


(Table  1  continued) 

Wenckebach  (Rg.  7c) 
pO  (SA  rhythm) 

pi  (AV  relative  ref.  sub-period) 
p2  (AV  relative  ref.  sub-period) 
p3  (AV  relative  ref.  sub-period) 
p4  (AV  conduction  delay) 
p5  (AV  conduction  delay) 
p6  (AV  conduction  delay) 
p7  (AV  conduction  delay) 
p8  (AV  absolute  refractory) 
p9  (AV  absolute  refractory) 
plO  (AV  absolute  refractory) 
pi  1  (AV  absolute  refractory) 
pi  2  (Vent,  absolute  refractory) 


pO  (SA  rhythm) 

pi  (Bundle  Kent  refractory) 

p2  (AV  refractory) 

p3  (AV  conduction  delay) 

p4  (Bundle  Kent  conduction  delay) 

p5  (Vent,  refractory  period) 

p6  (Vent,  refractory  period) 


Table  1 


(0.92,  0.96.1.00,1.04,1.08} 
{0.16} 

{0.16} 

{0.16} 

{0.16,0.20} 

{ 0.20,  0.24 } 

{ 0.24,  0.28 } 

{ 0.28,  0.32 } 

{ 0.48 } 

{  0.56 } 

{ 0.68 } 

{ 0.80 } 

{  0.28,  0.32,  0.36  } 


(0.92,  0.96,  1.00, 1.04,  1.08} 
{0.3} 

{0.3} 

(0.12,  0.14} 

(0.08,  0.10} 

{0.3} 

{0.3} 


